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1 Introduction 



In their programme of defining functions by nonlinear ordinary differential equa- 
tions (ODEs), Painleve and Gambier built a list |nj of fifty second order, first degree 
such equations, the general solution of which is uniformizable; among them, only 
six define new functions, called the Painleve transcendents (Pn). There have been 
recently many achievements |13| in obtaining some discrete analogues of these fifty 
equations. However, at present time, the nice features of the continuous case (de- 
gree one in u" for the ODEs, order two for the Lax pairs of the six (Pn) equations, 
coalescence cascade from (P6) to (PI) for the ODEs as well as for their Lax pairs) 
are far from being all implemented in the discrete case. Table || outlines some of the 
gaps to be filled. 



Table 1: Current state of the discrete (Pn) equations and their Lax pair. Column 

1 contains the type (difference or q— difference) of the discrete (Pn), column 2 the 
degree of this discrete equation in u, column 3 the order of the Lax pair if one is 
known, column 4 the appropriate references. A blank denotes a missing information. 
Ideally, there should be at least one entry for each d-(Pn), only digits 1 in column 

2 and digits 2 in column 3. 



d-(Pn) 
q - (Pn) 


degree 


order of Lax pair 


references 


d-(P6) 


2 




[El 


g-(P6) 


2 


2 


m 


d-(P5) 








g-(P5) 


1 




m 


d-(P4) 


1 




m 


d-(P3) 


2 


4 


p 




9 -(P3) 


1 


4 


m 2i] 


d-(P2) 


1 


2 




d-(Pl) 


I 


2 


i 




d-p 


I 


2 


m 





In this paper, we address the question of finding a second order Lax pair for each 
discrete (Pn) equation. 

In section [|, we present the five (Pn) equations chosen by Gamier, which have 
some advantages over the six usual ones. In section || we review the different 
kinds of continuous Lax pairs and discuss them according to their relevance for 
discretization. In section |], we convert the scalar "Lax" pairs given by Gamier 
in 1911 into an equivalent traceless matricial form, so as to remove the apparent 
singularity unavoidable in the scalar form. In section |5], we recall the discretization 
rules which we previously established. In section [], we give some preliminary results. 
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2 The five (Pn) equations of Garnier and their 
coalescence 

The (P6) equation for u(x) possesses in the plane of u four poles (oo, 0, 1, x) with 
the same residue equal to 1/2, and it depends on four parameters a,j3, 7,5. The 
scheme of the successive coalescences of these four poles 

(1/2,3/2) 



1/2,1/2,1/2,1/2) 



;i/2, 1,1/2) 



(2) 



(1,1) 

defines, from (P6)(u, x, a, j3, 7, 5), four other equations with four parameters, chosen 
as follows by Garnier |TT| 



(P6) u" 



1 



u u — 1 u — X 
u{u — l)(u — X 



u' 2 - 



1 



1 



x — 1 



(P5) u" 

(P4') u" 

(P3') u" 

(P2') u" 
(J') «" 



x 2 (x — l) 2 

1 1 

2u 4 

/2 



_ S X - 1 

a + + 7 7 — + 5 



U — X. 
x(x — 1) 



(u - ly 



(u — X) 



u' 2 -- + 

X 



u' 



(u - If 



u - 1 
u'" /3 o 

w' 2 m' aw 

3 , „„.\ , _yc.2 



era H — 
u 



u r u(u + V 
x u — 1 





Axu 2 + 2x 2 w) + 45(w 2 + xu) - 2au + - 
/ it 



4x 2 4x 4m 



5(2iT + rcu) + 7(6w 2 + x) + /3m + a 
2Su 3 + 67M 2 + (3u + a. 



Under the group of homographic transformations on u, the (P4') equation is equiv- 
alent to (P4) for 7^0; the (P3') equation is that defined by Painleve (Ref. p0| 
p. 1115) in the class of (P3), the (P2') belongs to the class of (P2) for 6 7^ and to 
that of (PI) for 5 = 0,7 7^ 0. This list has been terminated by the elliptic equation 
which is the autonomous limit of (P2'), which we denote (J') (like "Jacobi"). 

This is the most advantageous list of equations to be considered for the study 
of Lax pairs, for it allows to take full advantage of the four parameters a, /3, 7, 8 
appearing in each equation. 



The successive coalescences 



(P4') 



(P6) 



(P5) 



(P2') 



(J') 



(2) 



(P30 

from an equation E(x, u, a, j3, 7, 5) = to another equation E(X, U, A, B, C, D) — 
are described by the following Poincare perturbations 
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(x, u, a, (3, 7, 6) -> (X, U, A, B, C, D, e) 

6^5: (x, u, a, (3, 7, S) = (1 + eX, U, A, B, e~ x C - e~ 2 D, e~ 2 D), 
5^4': (x, u, a, (3, 7, 5) = (1 + eX, eU/2, 2Ce~ 4 + 28/^r 3 , 5/4, 

- ACe- 4 - Q0De~ 3 , 2Ae~ 2 - 2Ce~ 4 - 32De- 3 ), 
5^3': (x, u, a, (3, 7, 6) = (X, 1 + eU, e^A/A + s~ 2 C/8, -e~ 2 C/8, eB/4, e 2 D/8) 
4' -> 2' : (x, w, a, /3, 7, 5) = (e 2 X/4, 1 + el7, -4fi£- 4 + 96<^ 5 - 24L>£- 6 , 

lQAe- 3 - 8Be- 4 + A8Ce~ b - 8De~ Q , 

- 32Ce- Fj + lQDs- 6 ,48Cs- 5 - lQDe~ 6 ), 

3' -»• 2' : (x, m, a, /?, 7, 5) = (1 + e 2 X/2, 1 + el/, 16Bs~ 4 - QACe~ 5 - 32L>£- 6 , 

32£>£- 6 , -8Be~ 4 + 48Ce~ 5 + 16De~ 6 } 
lQAe' 3 - 8Bs- 4 + 16C^ 5 - 16L>£" 6 ), 

2' -»• J' : (x, u, a, (3, 7, 5) = (eX, e^U, e~ 3 A, e~ 2 B, e^C, D) 

in which e goes to zero. These six transformation laws on x and u are affine, contrary 
to those of the usual coalescence cascade. 

3 Various Lax pairs of (Pn) in the continuous case 

Lax pairs can be obtained either by the isomonodromic deformation of a given 
second order linear differential equation (defined either in scalar form or in matricial 
form), or by the reduction of the Lax pair for some partial differential equation 
(PDE) ||, which makes three possibilities. 

3.1 Lax pairs by scalar isomonodromy (Garnier) 

In order to find (P6), Richard Fuchs had to assume in the second order scalar 
equation 

d 2 ^ + {S/2)j> = 0, (3) 

the presence of four Fuchsian singularities, put a priori at the points t = 00, 0, 1, x. 
As shown by Poincare (Ref. |E§ pp. 217-20), the maximum number of apparent 



singularities (a singularity is apparent iff in its neighborhood the ratio of any two 
solutions i[) is singlevalued) which can be added in this case (second order, four 
Fuchsian singularities) is one. If no apparent singularity is added, a case which 



corresponds to the Heun equation [^6|, P7 |, the monodromy is trivial. By adding one 



apparent singularity, at a location denoted t = u, R. Fuchs || || found that u(x) 
had to satisfy the (P6) equation. This result has been completed by Garnier, who 
established a coalescence cascade starting from (P6) and yielding scalar pairs for 
(P5), (P4'), (P3'), (P2') and (J') while respecting a remarkable symmetry between 
t, x and u. 

These scalar Lax pairs of Garnier are linear in a, (3, 7, S and characterized by the 
two quantities (S,C), with the cross- derivative condition Z = 

d 2 ^ + (S/2)if> = 0, (4) 
d x ifj + Cd t iP - (1/2)CW = 0, (5) 
Z = S X + C ttt + CS t + 2C t S = 0. (6) 
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The value of S is given by 

S 3/4 /W + ft [{^u'f - (3 2 }e + f G (u) 
"2 = WW 2 + e 2 + /G(t) ' (7) 

in which the function /q and the various scalars are defined in Table ||. 



Table 2: Scalar Lax pairs of the Painleve equations 





(P6) 


fP5) 


\ L ^ ) 




\ L * ) 


Pi 


Tl T — 1 1 

Z(u — x j 


rv> 

z(u — 1) 


1/4 


rv> 
Jb 

ZU 


-1 




—u + 1/2 


— u + 1/2 


1/2 


-1/2 





fo(z) 


a b 

z 2 + {z-iy 

c d 


a 6x 

2 2 (^-l) 3 
cx 2 d 


a 8(z + 2x) 

z 2 4 
7(2 + 2x) 2 

16 


6x 

a 7 

cx 2 d 

+— + " 

z 


2az + (3z 2 
+2 7 (2z 3 + xz) 
+5{z i + xz 2 ) 


(z — x) 2 z(z — 1) 


' (2- l) 4 ' (z- l) 2 


a 


2(a + b + c + d+ 1) 


2(a + d+ 1) 


-46 


8d 






-2 (a + 1/4) 


-2(a+l/4) 


-8a - 2 


86 




7 


2(6+1/4) 


26 




16a 




5 


-2c 


-2c 




-16c 




e 


u — X 
u(u — 1) 


1/u 


1/it 


1 


1 


ei 


i(t-l) 


f (t - 1) 


-t 


t 


1 


e 2 


t(t - l)(t - x) 


t{t -l) 2 


t 


t 2 


1 


-C 


t(i - l)(u-x) 


t{t-l){u-l) 


t 

2 


tu 


1/2 


(t — u)x(x — 1) 


(t — u)x 


t — u 


(t — u)x 


t — u 



3.2 Lax pairs by matricial isomonodromy ( Jimbo and Miwa) 

Equivalently, the second order deformation equation can be replaced by a first 
order, two- component, matricial system. The advantage is the possibility to get rid 
of the apparent singularity, unavoidable in the scalar case, by a change of basis. The 
four Fuchsian singularities are similarly put at the a priori locations t = 00, 0, l,x. 



The resulting matricial Lax pairs obtained by Jimbo and Miwa |15| are different 
from those of Gamier once converted to matricial form; in particular, they are no 
more linear in cn,/3, 7, S, but algebraic. Some of them, namely those for (P3), (P5) 
and (P6), are restricted by nonvanishing conditions on the parameters. We give 
here traceless Lax pairs of (P3), (P5), (P6) for generic values of a, (3, 7, 8 (indeed, in 
Ref. some of these parameters are restricted to numerical values, and the Lax 
pair of (P6) is not traceless). In the following formulae, the symbols o~k denote the 
Pauli matrices 

° l = (l o) ,a2= { < i o0' a3= (o -1 ) ' aj(Tk = Sjk + ie i kl ° u ( 8 ) 
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a 



1 




a 





1 



while the symbol Z is the commutator of the Lax pair 

<),.c = L^, ihr = Mij), Z= <),. -L,dt~ M] 

and the symbol E denotes the (Pn) equation. 
The pair of (P3) @ for (7, 5) ^ (0, 0) is 



(dx/4 + t 2 (cx/4 - z))a 3 



M 



u 2 Z 



iz 
V 

w'/w 



+ (r 2 - r 1 u)wa + + (fV - t~ 2 z(z - cx/2))w~ x a' 

(dx/4 - r 2 (cx/4 -z)- r 1 ^ -1 - 2)/4)(7 3 

-(r 2 + r 1 u)wa + + (fV + t _2 z(,z - cx/2))w~ 1 a~ 

((l/8)r 2 (2xa 3 + (x 2 u- 2 (u' -d) + ((3d' 1 - l)xit -1 )w~V~) 

-(l/Sjr^^u -3 ^' - d) + x)w~ 1 (T-)E 

xu' /u 2 + cx — dx/u 2 + ((3d' 1 — l)/u, 

-z 2 u + z(cxu + (3d' 1 - 2)/2 + (a + 2c - Pcd~ x )x/%, 

-u'/u + d/u + (1 - (1/2)^0/^ 
7, d 2 = -5. 



The pair of (P5) 0] for 5 ^ is 

L = dta 3 /2 

+x~ l (z + O - u(z + (0 O - X + e oo )/2))wa + 
+x~ l (-z + u-\z + (0 O - 0i + e oo )/2))w"V 

?o + 600)72, 



M 



x z + 6o/2 z 

id — 

K 2 t 



t 



-)03 



+ ((^ + O )A - U(Z + (0 O - X + eoo)/2)/(t - 1))W 

+(-z/< + u~ 1 (z + (0 O + 0i + eoo)/2)/(t - i))«rV 

— it tt(t — 1 — tw)w \ X 

(t + u-tu)/w u ) t(t - l)u(u - 1) 

u'-du , (30 o + 0i + 600) - (0 O - 0i + Qoo)u 



E 



-X 



+ 



(9) 



(10) 

(11) 

(12) 
(13) 
(14) 
(15) 
(16) 



(17) 
(18) 



w'/w 
d 2 



(19) 
(20) 

2(u-l) 2 ' 4(u-l) ^ 

-2z - O + u(z + (0 O - 0i + 600)72) + vT\z + (0 O + 0i + 6oo)/2)]/x (22) 
25, (0 O - 0i + Boo) 2 = 8a, (0 O - X - 0oo) 2 = -8/3, (1 - O - 0i) = l/d. 



The traceless pair of (P6) for (2a — 1,2/3 + 1, 27 — 1,8)^ (0, 0, 0, 0) can be found 



in Ref. 18 



3.3 Lax pairs by reduction of a PDE 

One first selects a PDE which admits a second order Lax pair. The reductions 
create singularities at four points at most, which must then be moved by homography 
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to t = oo, 0, 1, x. The resulting Lax pairs are generically algebraic in a, (3, 7, 5, maybe 
with some restrictions on them. 

One such system is the pumped Maxwell-Bloch system, defined either in the 
complex form 

e T = P, Px = Ne, e T = p,p x = Ne, N x + (pe + pe)/2 - 4s = 0, (23) 

with s constant (the system is "pumped" when s is nonzero), or in the real form 
with one more arbitrary function m\(T) 



e T = P, p x = Ne + mi{T), N x + pe - 4s = 0. 



(24) 



One does not know a possible common parent to these two systems. These sys- 
tems describe phenomena in nonlinear optics, stimulated Raman scattering and 
self-induced transparency. Their respective Lax pairs are, for the complex system 



D [0 



e 
-e 



+ / 



for the real system 



e 
-e 



+ / 



1 
-1 



1 
-1 



1 



M 



1 

4? 



N -p 
-p -N 



(25) 



N 



M Af\-p-m l /{2f) 



-p + mi /(2f) 
-N 



(26) 



with j 2 = 2sT + A 2 and A an arbitrary complex constant. 

A Lax pair of (P3) is obtained by the following reduction of the pumped real 
Maxwell-Bloch system (M) E8L pi g 



x 
t 
L 
M 



XT 1/2 , e = -iT 1/2 [u'/u + cu + d/u], m x = -i{ad - (3c), s = -cd/2, (27) 
T- 1/2 f(T), (28) 
(l/2)(u'/u + cu + d/u)a 3 + ta u (29) 
x/t)L + t~ 2 cdxa 1 /2 + t~ 3 (-ad + /3c)a 3 /A (30) 

t 2 



r 2 (au + >yxu 2 + c(xu)')a + /2 + r 2 ((3/u - 5x/u 2 - d(x/u)')a'/2 

2 



Z 

„2 



±- { a,-r\cua^du^a-))^— d 
7 ,d 2 = -5. 



E. 



t 2 + cd } 
(31) 
(32) 



The transition matrix 



ft 1 ' 2 

I t- 1 / 2 



makes this pair invariant by parity on t. For 

7# = 0, it has two nonFuchsian singularities t 2 = 00, 0. For 7$ 7^ 0, it has three 
singularities in the variable —t 2 /(cd) : two Fuchsian and 1, one nonFuchsian at 
infinity, i.e. one more than expected for (P3), and there exists no coalescence of two 
of these three singularities which preserves (P3). Therefore this is in fact a pair of 
(P5) for 5 = |28| since (P5) for 5 = 0, 7 7^ is equivalent to the full (P3) under a 
birational transformation lO, 161 . 
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Another Lax pair of (P3) is obtained by a reduction of the unpumped complex 
Maxwell-Bloch system (^), but this one is restricted to (7,5) 7^ (0,0) 

x = XT 9 , e = T^- 20 ^,e = T^ 1+20 ^, s = 0, t = T- e f(T), (33) 
L = -(l/2)(u'/u+d/u+l/x)ai + (a/(cx) + cu)(a 3 + a 2 /i)/2-tux/c, (34) 
M = \t~ 2 cdiy~ l u' l (u/u + d/u-((3 + d)/(dx))+r 1 a/c}(a 3 + a 2 /i)/A (35) 

+ — vx 2 /ca^ + t' 1 (a/c)o'2/i — t~ x(u'/u + d/u + \/x)a\ — t~ 2 v (?d<j\ /4 

Z = (l/4)u -2 [t- 2 i/ -1 (<T3 + tr 2 /i) -t^xua^E, (36) 
^ = -t*/(4c) (37) 



-5. 



(38) 



in which is arbitrary. The arbitrary parameter z/ just reflects the scaling law of 
(P3) and it can be set for instance to c, d or 1. This pair is much simpler than the 
one of Jimbo and Miwa (|T0|)-( |TT|) , and it possesses two nonFuchsian singularities 
t = 00, 0. 

A Lax pair of the full (P5) is obtained by a reduction of the pumped complex 
system fl23|) [17] with the restriction 5^0 



x = XT 1 / 2 , e = T l ^F, e = T l ' 2+t ^G, s 
t = 2T~ l/2 f(T)/d+l/2 

(2t - l)4d F/2 
-G/2 -(2t-l)4d 

T " 



-5/16, 



(39) 
(40) 



L -. 
M 
Z 
F 



t- 1 

d 2 x 



{t-\){u-\f 

-1/2 



u 



G = u- 1 ' 2 



u 

u' + du 



(41) 

+ (d/2)xa 3 (42) 
[cr 3 + u~ l/2 {(t + u- tu)va + + (t - 1 - £u)v~V~)]£, 

(43) 



u' + du „ . , „ . . . 
— + 2(t c (it- 1) -tju)/x 



u- 1 

,2 1 i\ + /„,2 



— 2{t c {u — 1) — tju)/x 



-1 



= (^-(u 2 + 1) - i c (u 2 - l))/(x«) 
C + d = -t,a 3 - (ar/2) / ° j 



Cq — Ci 



-G 

a^ii' + du) (7 + d)(u+l) 



0-3 



P = 
d 2 



(u-1) 2 2d(u-l) 

+ n _1 ^ 2 ((p + q)va + + (p — q)v ~V~) 
x(m + !)(«' + tZtt) , (7 + d)u 



2(u-l) 2 



-u' + 



d(u- 1) 



q = t c {u + 1) — tjU 



-26, t 2 c 



-13/2, (t c - tj) 2 = a/2. 



(44) 

(45) 
(46) 

(47) 
(48) 
(49) 
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3.4 Comparison of these three kinds of Lax pairs 

Our goal is to discretize as easily as possible without introducing restrictions on 

The three kinds of Lax pairs are compared in Table |] according to several criteria. 

Table 3: Advantages and inconveniences of the three kinds of Lax pairs. The ideal 
situation would be that all entries be "yes" . 





scalar Gamier 


matricial Jimbo and Miwa 


by reduction 


no apparent singularity 


no 


yes 


yes 


linearity in a, /3, 7, S 


yes 


no 


no 


easy confluence 


yes 


no 


yes 


arbitrary a, (3, 7, 5 


yes 


no 


no 



We therefore choose the scalar pairs of Gamier and, before attempting to dis- 
cretize, we remove the apparent singularity by converting them to matricial form. 

Remark. If one first discretizes a PDE then performs a discrete reduction, the 
resulting discrete ODE is not explicit but is defined by the system made of the 
discrete PDE and the discrete nonlinear relation defining the reduction; for details, 
see section 2.4 in Ref. fli"3|. 



4 The Garnier matricial Lax pairs 

From the scalar Garnier pairs ©-(0), let us build an equivalent matricial pair 
which does not contain any more the apparent singularity t = u. 
Let us choose the two matrices as traceless 

-t))(*H (50 » 

A tf))(t)^. (51) 



VM 21 -MuJJKip: 
The elimination of one component, say ipi, defines the scalar system 

d 2 t - + M n ,t ~ M 2 U - Mn^ - M 12 M 21 ) ^ = 0, (52) 

d x - ^d t + L u - M u ^f) ^ 2 = 0. (53) 
M21 M21 / 

The coefficients of this scalar Lax pair have simple poles at M 2 \ = 0. Let us first 
convert the double pole t = u of the coefficients of the pair of Garnier into a simple 
pole, by the transformation 

$ = y/t - ui/j, (54) 

( s ' 2 -^ + (f + (^))^°- < 55 » 
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with S given by (0) and C by Table 0. The six coefficients (Ljk, Mjk) of the Lax 
pair are obtained by identifying the singularities M 2 i = and t — u = in (|52"D-(|5"3|) 
and (E5|)-(E3). This results in 



M u = 


/W + 0() 

J 




M 12 = 


M 1M - Mi! - 


(t - u) 2 


M 21 = 


t — u, 




L21 = 


-(t-u)C, 




L n = 


u' - 

CM n + 


(C + (t- 
2(t - u) 


L\2 = 


-L lltt + M 1M 


+ L 2 iM 12 


t — u 





'S 3/4 



t-u\2 (t-uf 



(57) 

(58) 

(59) 
(60) 

(61) 
(62) 

Once expressed in terms of (t, x, u, u', a, 0, 7, 5) only, this matricial Lax pair has 
the following features : regularity at t — u, concentration of the dependence on 
a,0,j,5 in the elements (L 12 ,M 12 ), linearity of these two elements in a,0,j,5, 
absence of any restriction on a,0,j,5. For instance, the element Ln evaluates to 

£11= ,0,0,0 (63) 

(x(x — 1)) X 

for the five successive (Pn) of Gamier, and the element L 12 to 

Ln = - TT , 0, 0, 0, (64) 

t{t — l)(t — x)u(u — 1) 

for a = b = c = d = 0. The confluence operates on each element separately. 
For instance, the matricial Lax pair of (P2') is 

L =(l/2 o) + (o J)[/ 5 + 4 7(^ + 2w)+5(t 2 + 2t M + 3n 2 + a;)]/2, (65) 
M=(~* °/) + (° l)[2a + 0(t + u) +1 (4(t 2 + tu + u 2 )+2x) 



+ S(t 3 + t 2 u + tu 2 + u 3 + x(t + u))), (66) 



and the one of (J') is 



L= (l/2 o) + (o J)[^ + 4 7(A + 2n) + 5(A 2 + 2An + 3n 2 )]/2, (67) 
M=(~ U ' °,) + (° n l)[2a + 0(X + u) + A 1 (X 2 + Xu + u 2 ) 



X — u u'J \0 

+ 5{X 3 + X 2 u + Xu 2 + u 3 )}, (68) 

in which A is an arbitrary constant. 
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5 Rules of discretization of matricial Lax pairs 



In order to discretize a given ODE having the Painleve property, some basic 
discretization rules have been proposed 

The definition of the (continuous) matricial Lax pair (L, M) can be discretized 

as 

ip(x + h/2) = A{x, z)tl)[x - h/2), d z ip(x - h/2) = B(x - h/2, z)ip(x - h/2), (69) 
d z A(x, z) + A(x, z)B{x - h/2, z) - B(x + h/2, z)A{x, z) = 0. (70) 

The continuum limit /i — ^ is then 
A - I 



L, (dz/dt)B^M, 



h _ 

(dz/dt)(d z A + AB- BA)/h -»■ d t L - d x M + LM - ML, (71) 

with some link F(t, z,h) = between the spectral parameters t and z. For a 
second order equation E(u,u,u,x,h) = 0, the operators A and B must have the 
u— dependences A(u, u, u), B(u, u). 

Recently, we proposed || a direct method for finding the Lax pair of a discrete 
equation whose continuous pair is known. This consists in discretizing the continuous 
Lax pair by obeying the following common sense rules 

1. conserve the matricial order. This is indeed the differential order of the scalar 
Lax pair, which must be conserved; 

2. replace the continuous spectral parameter t by an unspecified function T(z, h); 

3. discretize the operator L centered at the three points x — h, x, x + h. If L is 
traceless, so is its discretization; 

4. discretize the operator M centered at the two points x — h,x. If M is traceless, 
so is its discretization; 

5. replace each monomial (du/dx) k by its discretization obeying the general rules, 
multiplied by the fc-th power of an unspecified function g(z, h). This function 
g, whose continuum limit must be 1 for any z, represents the ratio of the 
stepsize h to the differential element dx; 

6. take B as the product of the discretized M by an unspecified function J(z, h) 
(like Jacobian) representing a discretization of the derivative dT/dz; 

7. take (A — l)/h equal to the sum of the discretized operator L and a diago- 
nal matrix of unspecified functions of (z, h) only, diag(gi, #2); these functions, 
whose continuum limit must be zero, account for the dissymmetry of the for- 
mula defining A. 
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6 Application to the discretization of the (Pn) 
equations and of their Lax pair 

The (Pn) equations and the pairs of Gamier, in their scalar form as well as 
in their matricial one, are linear in the four parameters a,/3, 7, 5. If one requires 
this property to be conserved in the discrete case, the problem of finding a discrete 
equation and its discrete Lax pair splits into five successive subproblems : all four 
parameters zero (or numerical constants), coefficient of a alone, of /3 alone, etc. 

When applied to the autonomous limit (J') of (P2') and its Lax pair (|67| )-(l68D , 
the discretization rules give straightforwardly the desired result, namely 

u — 2u + u , r 2 ,_ 



+ 5u (u + u) + 2^u(u + u + u) 
h l 



+ /?(A 1 « + (l-A 1 )^)+a = 0, (72) 



F(x + V2)-F(x-V2) 
u(x + h) - u{x - h) 

B/J=(^~- ) (i , ° w*,Wn l)[2a + (3(T+(u + u)/2) 
' \T - (u + u)/2 g(u-u)/hj \0 0J l HK K ~" ' 

+ 7 (4T 2 + 2T(u + u)+ Auu) + 5F(x - h/2)] , (74) 

_ T 4 - (u(x + h/2)u(x - h/2)) 2 - T 2 {u{x + h/2) - u(x - h/2)) 2 /2 

(X} ~ T - {u{x + h/2) + u{x - h/2))/2 ^ ' 

s(^zi = ^ +7r+ 4 (76) 

T'(z) = 0. (77) 

In these expressions, Ai is an arbitrary constant (which could be absorbed in the 
discrete second derivative), the functions J(z), T(z) and g(z) are constant, the 
spectral parameter is the constant value of T. 

One notices an unexpected drawback : the apparent singularity t = u, absent 
from the continuous matricial pair, reappears in the coefficient of 5, because the 
polynomial T 3 + T 2 u + Tu 2 + u 3 is discretized as the quotient of a discretization of 
T 4 — it 4 by a discretization of T — u. 

This is the reason why the Lax pair of d-(P2') obtained by this method is, for 
the moment, restricted to 5 = i.e. to d-(Pl) 

u — 2u + u „, 2 



+ 8{u\u + u) + xu) + 7(2u(m + u + u) + x) + f3u + a = 0, (78) 



h 2 

^=(1/2 9 %) + (°o 1)W + W + *)], P9) 



T-(u + u)/2 g(u-u)/hj \0 0. 

+ 7 (4T 2 + 2T(u + u)+ Auu + 2{x - h/2))], (80) 

9{Z) 1~ 1 = /3/4 + 7 T, (81) 
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1 T\z)= 1 J{z)g{z), 
g\z) = (l/2)jh 2 J(z), 



(82) 
(83) 



with J(z) arbitrary. The contributions 5G(x 
A12 and Byi lead to an impossibility. 

For (P3'), the result given in Q deals with the case a 
the degree of the d-(P3) given in Ref, 



and 8F(x — h/2, z) in the elements 



this case a second order discrete Lax pair can be obtained. 



P = 7 = 5 = where 
12] reduces from two to one (see Table 



; m 



7 Conclusion 

This direct approach, which serializes the difficulties, can be improved in the 
following directions. If one sticks to the matricial form, one should design more 
efficient discretization rules, so as to at least find the contribution of 5 in the Lax 
pair of d-(P2'). The other possibility is to directly discretize the Gamier pairs on 
their scalar form fl55|)-(|56D . 
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